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ABSTRACT 

We consider grains with superparamagnetic inclusions and report two new condensed matter effects 
that can enhance the internal relaxation of the energy of a wobbling grain, namely, superparamag- 
netic Barnett relaxation, as well as, an increase of frequencies for which nuclear relaxation becomes 
important. This findings extends the range of grain sizes for which grains are thermally trapped, 
i.e. rotate thermally, in spite of the presence of uncompensated pinwheel torques. In addition, we 
show that the alignment of dust grains by radiative torques gets modified for superparamagnetic 
grains, with grains obtaining perfect alignment with respect to magnetic fields as soon as the grain 
gaseous randomization time gets larger than that of paramagnetic relaxation. The same conclusion 
is valid for the mechanical alignment of helical grains. If observations confirm that the degrees of 
alignment are higher than radiative torques can produce alone, this may be a proof of the presence of 
superparamagentic inclusions. 

Subject headings: polarization -dust extinction -ISM: magnetic fields 



1. INTRODUCTION 

Alignment of dust in magnetic fields is one of the as- 
trophysical problems of the longest standing. Discovered 
by Hiltner (1949) and Hall (1949) via observation of the 
polarization of starlight attenuated by interstellar dust, 
aligned grains are shown to trace magnetic fields in vari- 
ous environments, from interplanetary space (see Rosen- 
bush et al. 2007) to circumstellar regions (see Aitken et 
al. 2002; Cho & Lazarian 2005). 

The presently-accepted dominant alignment process is 
based on the action of radiative torques (RATs) on irreg- 
ular grains (see Lazarian 2007 for a review). This process 
was proposed by Dolginov & Mytrophanov (1976) and 
extended, corrected and elaborated in the subsequent 
publications (Lazarian 1995; Draine 1996; Draine & 
Weingartner 1996, 1997, henceforth, DW96 and DW97; 
Weingartner & Draine 2003, Weingartner 2006; Lazarain 
& Hoang 2007a, henceforth LH07; Hoang & Lazarian 
2008). This has induced theory-based modeling of po- 
larized radiation from molecular clouds (Cho & Lazarian 
2005; Pelkonen et al. 2007, Bethell et al. 2007) and 
accretion disks (Cho & Lazarian 2007). 

While DW96 professed the omnipotence of RATs for 
aligning interstellar grains, the physics involved happens 
to be more complex than it was initially anticipated. In 
particular, an analytical model (AMO) for RATs pro- 
posed in LH07 shows that, for a variety of the directions 
between the magnetic field and the radiation flux, the 
alignment is not perfect (see also a particular example of 
this in Weingartner & Draine 2003 and more details in 
Hoang & Lazarian 2008). This causes the justified wor- 
ries whether all the essential physical processes has been 
accounted for to explain the high degrees of interstellar 
polarization (see Whittet et al. 2008). 

In LH07a we showed that ordinary paramagnetic re- 
laxation does not affect the RAT alignment in most as- 
trophysical environments, e.g., in interstellar gas. This 
does not preclude inclusions with strong magnetic re- 
sponse, e.g. superparamagnetic inclusions, to affect the 



RAT alignment. Such grains were discussed in Jones & 
Spitzer (1967) and also Martin (1995). Given the abun- 
dance of iron in dust, this would not be a far fetched sug- 
gestion. Moreover studies of interstellar grains caught in 
the Earth atmosphere (see Bradley 1994; Goodman & 
Whittet 1995), as well as, laboratory studies of the sili- 
cate films (see Djouadi et al. 2007) are suggestive of the 
abundance of such inclusions. Below in §2 we consider 
both the effect of superparamagnetic inclusions on grain 
internal relaxation, i.e. on the relaxation that may in- 
duce grain to rotate about its axes of maximal inertia 
(Purcell 1979; Lazarian 1994; Lazarian & Roberge 1997; 
Lazarian & Efroimsky 1998; Lazarian & Draine 1999ab, 
henceforth LD99ab) and on the RAT alignment (in §3). 

2. BARNETT AND NUCLEAR RELAXATION IN 
SUPERPARAMAGNETIC GRAINS 

A general susceptibility of the grain with superparam- 
agnetic inclusions is given in Draine & Lazarian (1999). 
If, however, we consider only low frequency response 
of the material, superparamagnetic inclusions containing 
N c i atoms act in the external magnetic field as magnetic 
moments of the moment N c i/j,, where [i is the moment 
of an individual atom. This increases the grain magnetic 
susceptibility by a factor of N c i, which may vary from 
approximately 10 to 10 6 (see Kneller & Laborsky 1963; 
Billas et al. 1994). The enhanced susceptibility affects 
the internal relaxation through the Barnett effect. 

Barnett effect (Barnett 1915) can be well understood 
classically (see LD99a for a quantum mechanical treat- 
ment). If one considers, instead of an electron spin, a 
gyroscope with a momentum J. As the grain rotates 
with angular velocity f2, the torque acting on a spin is 
2jJ = O x J. The equivalent torque can be induced 
by a magnetic field Hg acting along the grain angular 
velocity axis, i.e. 1/c/ixHs, which provides the Barnett- 
equivalent field H eqv — ^f2. The magneto-mechanical 
ratio for an electron is 7 = /i/( Jc) = e/mc, where e and 
m are the electron charge and mass, respectively. Thus 
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the Barnett-equivalent field is H eqv = 

The Barnett-equivalent magnetic field induces grain 
magnetization M = x^-eqv = X^/l- Purcell (1979) no- 
ticed that when a grain wobbles, the precession of J~2 in 
grain axes induces the precession of the vector of magne- 
tization M, which causes energy dissipation. For a grain 
at an absolute zero temperature this would result in the 
grain getting to rotate about its axis corresponding to the 
maximal moment of inertia. For grains at finite tempera- 
tures, the relaxation induces important coupling between 
rotational and vibrational degrees of freedom, which re- 
sults in many important effects (Lazarian 1994, Lazarian 
& Roberge 1997, LD99a,b). We expect for grains with 
enhanced magnetic susceptibilities, i.e. larger \ both the 
magnetization and the relaxation that this magnetization 
entails to be enhanced. 

To quantify the effect, consider the internal relaxation 
associated with the Barnett magnetization of a wobbling 
oblate grain. The angle 9 between the axis of major in- 
ertia ai and the angular momentum J changes according 
to LD99a as 

^ = -(h- 1)^-^( X " Mshx0cosdf(h,9), (1) 

where 7|| and I±_ are the moments of inertial parallel and 
perpendicular to the axis of major inertia, h = Iu/I±, 
oj = (h — l)Jcos0/7|| is the precession frequency, V is 
the grain volume, x is the imaginary part of magnetic 
susceptibility, and f(h, 6) is a function of order unity 
that takes into account details of spin-lattice coupling. 
For sufficiently slow rotation, e.g. < 10 6 this function 
is identically unity and Eq. (fl]) coincides with the result 
in Purcell (1979). 

If the magnetic moment of electrons on the outer 
partially- filled shell is p/j-b, where \lq = eh/2m e c is the 
Bohr magneton, while for bulk Fe we takep = 2.22 (Mor- 
ish 1980), then for zero-frequency susceptibility of a grain 
is (see Draine 1996) 

Xsu P (0) « 1.2 x lO- 2 N cl f p h 23 plf{- 5 \ (2) 
where f p » 0.01 is the fraction of Fe atoms, 77-23 = 

/(10 23 cm~ 3 ) and T15 = T/(15 K) are normalized grain 
density and temperature, respectfully. Superparamag- 
netic clusters undergo thermally activated remagnetiza- 
tion at a rate 

t-J « vo exp [-N cl T act /T] (3) 

where v Q « 10 9 and T act » 0.01LK" (see Morish 1980) 

Adopting critically damped susceptibility from Draine 
& Lazarian (1999) one can write the imaginary part of 
the grain susceptibility as x" = Xsup(0)wr / (l + (wr /2) 2 ) 2 
. Combining the latter expression for x with Eqs. (JTJ) 
and ([2]) one obtains for an a x a\/3 x a\/3 brick the time 
of alignment equal to 

TB „„,»ioV^(f) 2 x(i + (^) 3 ) 2 , 

(4) 

which is a factor ~ N c i smaller than the usual Barnett 
relaxation time. 

The competing process of internal relaxation, namely, 
"nuclear relaxation" was introduced in LD99b. While 



the magnetization arising from nuclear spins of a rotat- 
ing grain is ~ 1000 times smaller than the magnetiza- 
tion arising from electron spins, the nuclear relaxation 
~ H 2 qv is approximately ~ 10 6 times faster than the Bar- 
nett relaxation in normal paramagnetic material. This 
difference in the value of Barnett-equivalent field H eqv 
can be understood if we recall that the flipping of both 
electron and nuclear spins in order to get spins mostly 
parallel to grain angular velocity is entirely mechanical 
effect. For a given temperature the equilibrium excess 
of oriented spins depends only the angular momentum 
of the spin and the temperature (see Lazarian & Draine 
1999a). For both nuclear and electron spins the former 
values are equal of h. However, to get the same excess 
of orientation acting by external magnetic field on much 
weaker nuclear moments, one has to use a stronger mag- 
netic field. 

Superparamagnetism does not change the amplitude 
of the nuclear relaxation rate, but, nevertheless, affects 
the rate of the wobbling at which the relaxation gets 
saturated. Indeed, the relaxation of the system of spins 
depends on the rate of spin precession in the internal 
varying magnetic field of the neighboring magnetic mo- 
ments. In the medium with magnetic moments \x m with 
density n m one can use van Vleck (1937) expression 
for the internal magnetic field Hi ss 3.8n m /i m . If the 
medium is extended, Hi does not depend on whether the 
magnetic moments are clustered or not, as n m ~ A^ 1 , 
while (j, m ~ N c i- If T~ p is much faster than the time 
of nuclear moments precession uin gN^^Hi/h, where 
/Ltjv = eh/rripC is nuclear magneton, then the fluctuating 
magnetic field will induce a random walk with the step 
~ T sp . The change in direction of the nuclear moment 
per step is SO ~ ujnt sp , which means that the initial 
phase relationship is lost after a time 

Tnsp ~ (l/S6) 2 T sp « Ujjj 1 /(w N T sp ), (5) 

which should be compared with electron-nuclei decorre- 
lation time r„ e w Lu~^ l (Lu e /u)N) and a comparable rate of 
nuclei-nuclei decorrelation t„„ « 0.58T ne (77 e /?7„), where 
n e and n n are densities of unpaired electron and nuclear 
spins, respectively (LD99b). Therefore, if the rate of the 
nuclear relaxation T~Jp can be increased in the presence 
of superparamagnetic inclusions by a factor r sp Lu e , which 
scales according to eq. ((3} as exp (N c iT act /T) uj e , pro- 
vided that, nevertheless r~ p l 3> ujn- This means that 
for clusters with N c i < 1.5 x 10 4 the extension of the 
range of frequencies for the dominant nuclear relaxation 
exist. If clusters are larger, then superparamagnetically 
enhanced nuclear relaxation gets suppresses by the factor 
(1 + (ujT nsp /2) 2 ) 2 similar to eq. J4]). Figure [1] illustrates 
the relative roles of the Barnett and nuclear relaxation 
processes in paramagnetic and superparamagnetic grains 
and shows that for N c i = 10 4 , nuclear relaxation get 
modified in the range of grain sizes smaller than 2 x 10~ 5 
cm. As a result, for grains in the range from 3 x 10~ 6 
to 2 x 10~ 5 cm superparamagnetic Barnett relaxation 
dominates. 

The actual relaxation rate is important for getting the 
rates of grain flipping. The corresponding expression ac- 
cording to LD99a is = exp [,J 2 j 3 2 d — 0.5) where t~J 
is the internal relaxation rate, which is the sum of the in- 
ternal relaxation rates acting on the grains, J is the mag- 
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Fig. 1. — Time scales of internal relaxations: Barnett relaxation 
(Tflor), Barnett superparamagnetic relaxation (TB ar sup ), nuclear 
relaxation (r nuc i) and internal relaxation time (r;„ t ) as functions of 
grain size, for two values of angular momentum J/Jd = 5 {upper) 
and J/Jd = 1 (lower). Here we adopt number of iron cluster 
N c i = 10 4 . Relaxation by superparamagnetic inclusion is dominant 
for grains smaller than 3 X 10 cm for the case J/Jd = 1. 
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Fig. 2. — Trajectory maps for the alignment by RATs only 
(upper panel) and by RATs plus superparamagnetic relaxation with 
tDG,sup/tgas = 1 (lower panel) for Shape 4* and one light direction 
ifi = 30°. Unlike in LH07, only trajectories corresponding to J 
parallel to grain axis a are shown. Repellor points A in the upper 
panel becomes the attractor point A' in the lower panel as a result 
of superparamagnetic alignment. 



nitude of the angular momentum, and Jd — v^AsTd 
is the angular momentum of dust grain corresponding to 
the dust temperature T^; I\ is the moment of the axis of 
major inertia ai and ks is the Boltzmann constant. 

3. RAT ALIGNMENT OF SUPERPARAMAGNETIC GRAINS 

As we see in Figure [5] the DG relaxation within su- 
perparamagnetic material can influence grain dynamics. 
The maps shown follow the evolution of an irregular grain 
(Shape 4* in LD07, i.e., modified from Shape 4 by chang- 
ing the helicity from left to right), keeping track of its 
angular momentum J and the angle £ between J and 
the magnetic field. For simplicity, we disregard thermal 
wobbling and flipping (see Hoang & Lazarian 2008 for 
more details). Comparing the panels of Figure^ we see 
that without the DG relaxation (or in its presence, but 



assuming that grains are paramagnetic), only lower at- 
tractor point B is being observed in the phase trajectory 
map, while stationary point A is a repellor point (upper 
panel). However, this repellor point A becomes a new 
high attractor point A' when the grain is superparamag- 
netic (lower panel). 

When gaseous bombardment is included, then, as 
shown in Hoang & Lazarian (2008), grains are being 
kicked out from the low-J attractor point and end up 
at the high- J attractor point. Other randomization pro- 
cesses, e.g. shaking of a flipping grain by B x v elec- 
tric field arising from grain drift (Weingartner 2006) or 
from grain emission, plasma-grain interactions etc. (see 
Draine & Lazarian 1998) will have a similar effect. As a 
result, grain alignment with R w 1, i.e. perfect alignment 
of grains in respect to magnetic field, is achievable. 

What are the degrees of the paramagnetic enhance- 
ment that are required to create high- J attractor points? 
The answer depends on the angle between radiation di- 
rection and the magnetic field if) and grain properties, 
that, according to LH07a, can be conveniently parame- 
terized by the ratio q max = Q™ x ax / Q™ 2 al \ where Q™ x ax is 
the amplitude of the component of the torques parallel 
to vector k and Q"2 ax is the amplitude of the component 
of the torques in the plane of grain axis of major inertia 
ai and k. Let us estimate roughly when there are high- J 
attractor points in the presence of superparamagnetic in- 
clusion. For simplicity, we consider the case ■0 = 0°, and 
the default AMO with the mirror at an angle a = 7r/4. 
Thus the aligning torque reads (see eq. 25 in DW97) 
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where the first term containing a factor M = 
Urad 2 j°^ ftBa3 , where u ra d is the radiation energy den- 
sity, A is the wavelength, a e ff the effective grain size 
defined as the radius of the spherical grain of the same 
volume, t gas is the gaseous damping time, denotes the 
aligning component of RATs (see eq. 38 in LH07) and 
J' = J jI\u)T, where I± is the maximal moment of iner- 
tia and lot is thermal rotational rate at T — T gas . The 
second term represents the contribution from the super- 
paramagnetic relaxation with characteristic time 
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It is obvious from Eq. (J5|) that the position of the 
stationary points £ s = and tt (i.e. (F)^ = 0) is not af- 
fected by the presence of the additional relaxation. Also, 
the magnitude of spin-up torque (H)^ at this stationary 
point is not affected by the superparamagnetic relaxation 
either (see eq. 26 in DW97). As a result, the stationary 
point £ s is an attractor point when 



1 d{F)^ 



H) 



< 0. 



(8) 



At the stationary point £ s = 0, (H)$ > 0, therefore, 
using Eq. (J6J) and ||5J) one gets 
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DG.sp M (-Qf^ + 2Q™ 2 ax ) {-Q'l\ ax + 2Q™ ax y' 

where we have used J s /M = = O)) = Q^ ax 

is the spin-up torque at the stationary point £ s = 0. 
For Qff > 2Qf 2 ax , the criterion $LU§ is always ful- 
filled, i.e. there is high-J attractor points regardless 
of DG relaxation. This is similar to what obtained in 
LH07. When Q^ ax < 2Qf 2 ax , Eq. (E) gives t' DGsup < 
q max /(2 - q max ) with q max = Q™f x / Q™ ax . Thus, the 
AMO predicts that for q rnax — 0.8, the enhancement of 
paramagnetic susceptibility necessary to convert the re- 
pellor points into attractor points is t' DG sup < 0.67 for 
ip = . The same value q max corresponds to the Shape 
4* grain at A = 1.2/j.m and for this grains we get the nec- 
essary value that we obtain numerically t' DG sup « 0.56, 
which is close to what the AMO gives us. For experi- 
mentally studied in LH07 grain shapes the correspond- 
ing q max varied from 0.4 to 10 and for ip in the range 
from 0° to ~ 90° and found that the minimal value 
^'dc, sup ~ 5 x 10~ 2 is sufficient to produce the high- 
J attractor points in all cases. This is easily achiev- 
able for grains having superparamagnetic inclusions with 
N c i ~ 5 x 10 3 iron cluster. It is clear that such superpara- 
magnetic inclusions can substantially enhance the degree 
of grain alignment. 

To distinguish the RAT alignment in which magnetic 
field is only important in terms of inducing fast Larmor 
precession, from a hybrid type alignment that includes 
an additional process, namely, superparamagnetic relax- 
ation we shall call the latter process superamagnetic RAT 
alignment, or, in short, SRAT alignment. Note that su- 
perparamagnetic inclusions also enhance the rate of Lar- 
mor precession of J in respect to B, which makes the 
alignment with respect to B preferable to the radiation 
vector k (see discussion in LH07a). This may allow more 
reliable interpretation of polarimetry in terms of mag- 
netic field in circumstellar regions. 

4. DISCUSSION AND SUMMARY 

In terms of the RAT alignment it is possible to show 
that even inclusions with TV ~ 10 3 can result in the ap- 



pearance of a new high-J attractor point. As a result, 
the RAT alignment can be enhanced substantially, be- 
coming nearly perfect. This alignment is different from 
the alignment described in Draine & Weingartner (1996), 
however. In that paper, it was assumed that RATs pro- 
vide the spin-up similar to pinwheel torques in Purcell 
(1979), while ordinary paramagnetic dissipation does the 
job of alignment. The latter interpretation of the RAT 
action is erroneous, as RATs provide the alignment of 
their own, making the weak paramagnetic torques ir- 
relevant. It is only after being enhanced through the 
presence of superparamagnetic inclusions, that the para- 
magnetic dissipation can affect the alignment. If obser- 
vations prove that perfect alignment of grains is required, 
SRAT alignment can achieve this. Interestingly enough, 
the synthesis of pinwheel torques and superparamagnetic 
inclusions proposed in Spitzer & McGlynn (1979), does 
not work for interstellar medium, as superparamagnetic 
grains are strongly thermally trapped and cannot rotate 
suprathermally, unless RATs are present. Note, that 
superparamagnetic inclusions will also enhance the me- 
chanical alignment of helical grains. Such an alignment 
was proposed in Lazarian & Hoang (2007b) where it was 
shown that torques arising from the gas-grain drift may 
act on irregular grains the same way as RATs do. 
Our results above can be briefly summarized as follows: 

• Superparamagnetic inclusions affect both Barnett 
and nuclear relaxation, but in a different way: Barnett 
relaxation is being enhanced by the increase of mag- 
netic susceptibility, while the nuclear relaxation does not 
changes its amplitude, getting, nevertheless extended for 
higher frequencies. 

• Supraparamagnetic inclusions modify the RAT align- 
ment and RAT-superparamagnetic mechanism, which we 
termed SRAT alignment, may provide higher degree of 
alignment compared to the RAT mechanism alone. 

• The results above are also applicable to mechanical 
alignment of irregular grains. 
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